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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
TECHNICAL MEMORANDUM 1290 

THE STRUCTURE OF AIRY'S STRESS FUNCTION IN MULTIPLY 

CONNECTED REGIONS* . 

By Giusippe Grioli 


The determination of the stress state in a plane, homogeneous, and 
isotropic elastic system with assigned forces on the boundary and in 
the absence of body forces^ is known to lead, for the cases of both 
plane forces and strains, to the investigation of a biharmonic function, 
namely, the Airy function. 

With the aid of this function, various investigations have been 
conducted, particularly on elasticity problems of plane simply or doubly 
connected stress systems. 

In attempting to establish a systematic treatment of the problem 
of elastic equilibrium of a plane homogeneous and isotropic system of 
any shape and possessing any order of connectivity, the first step is 
to determine the singular terms (absent in a simply connected system) 
of Airy function in the most general case of stress and strain; on the 
basis of their significance, it may be decided which of these terms 
remain in various applied cases. 

It is therefore proposed to isolate the singularities of the Airy 
function for a general plane system and to show how these functions are 
connected by their mechanical properties with the state of strain and 
stress of the elastic system. 


*"Struttura della funzione di Airy nei sistemi molteplicemente 
connessi." Giornale di Matematiche, Vol. 77, 1947, pp. 119-144. 

!ln cases of interest in which body forces are present, as, for 
example, for constant body forces (gravity forces), the problem maybe 
reduced to one without these forces by changing the unknown function. 
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A special decomposition of the Airy function is therefore given; 
for simplicity of presentation, the following two different cases will 
be distinguished: 

(a) systems subjected to a plane dislocation and therefore to a 
state of stress with zero forces on the boundary 

(b) systems subjected to an external stress of a general type, but 
in the absence of dislocation (one-valued displacements) 

It is evident that from the two compositions corresponding to cases 
(a) and (b) the one is immediately obtained that corresponds to the 
general case of a plane system in which the state of deformation is the 
resultant of those deformations due to a general external stress and a 
plane dislocation. 

For a plane dislocation, only the components of the displacement 
parallel to the plane of the system are considered herein. Such a dis- 
placement results when, after the cuts that decrease the order of con- 
nectivity of the system are made, one of the faces of each single cut 
is given a displacement relative to the other faces. The displacement 
is composed of a translation parallel to the plane of the system and a 
rotation about an axis perpendicular to this plane; a certain quantity 
of material is then added or subtracted, as the case may be^. 

Treatment of the problems of plane dislocations with the aid of the 
Airy function seems natural, because the determination of the states of 
strain and of stress of a system subjected only to dislocation leads 
essentially to the problem of the integration of the equations of elas- 
tic equilibrium with assigned forces at the boundary. 

In the two decompositions presented, some of the singular parts 
figure directly as known elements in that their coefficients may be 
expressed in terms of assigned elements characteristic of the state of 
stress and strain to which the system is subjected. That is, among the 
coefficients of the singular terms, the characteristics of the dis- 
location figure in case (a); whereas in case (b), the components of the 
resultant forces and the resultant moments of the external stresses 
acting on the boundaries of the individual holes are involved. 


^Results on the plane dislocations of a circular ring will be found 
in references 1 and 2. 
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For the first case, the expression given by the Airy function satis- 
fies the condition of many valuedness of the required displacements from 
the dislocation considered; the condition of single -valuedness of the 
displacements is satisfied in the second case. The analytical problem 
that remains to he solved therefore reduces in each case to a boundary 
problem similar (but not identical) to that to which simply connected 
systems are reduced. 

It was also desired to show that even those coefficients of the 
singular terms that are not directly known are expressible in terms of 
characteristic elements of the deformation and that their invariance 
with respect to certain curves, with reference to which they are defined, 
corresponds to mechanical properties common to all continuous systems 
at rest, the known properties of the mean stress. 

In order to determine the character of the single- valued part of 
the Airy function in the two cases considered, a special decomposition 
of the biharmonic functions is employed in the manner of Poincare. This 
decomposition was established several years ago by G. Fichera (refer- 
ence 3) and was found to be very useful for the purpose. 

The many- valued part was determined by direct considerations, but 
it is necessarily present if, and only if, the acting forces on each 
singular boundary do not- constitute a system in equilibrium. 

From the given decomposition for the Airy function, those components 
of the displacement and of the rotation from which the singularities 
appear are derived without difficulty. It is thus possible to treat 
systematically also problems concerning systems in which the displace- 
ments as well as the forces are assigned on the boundary or partly the 
displacements and partly the forces (mixed problems) . Such decomposition 
formulas are given subsequently. 

For simplicity, a doubly connected system will first be presented; 
the results will then be extended to a system of any order of con- 
nectivity. 

The essential concepts of the theory of elastic dislocations that 
is used herein are directly derived from the fundamental investigations 
of Vito Volterra^. 


^Published in various notes in the reports of the R. Accademia dei 
Line ei in the years 1905-1906 and collected in one memoir in reference 4. 
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The author wishes to express his sincere thanks to Signora Virginia 
Volt err a, widow of Vito Volterra and to his son Prof. Enrico Volterra 
who kindly furnished him with the Italian and French editions of the 
papers . 


1. Many-Valuedness of Airy Function 

Assume the plane system defined in a region T hounded externally by 
the curve Co* internally by the curve Cp and referred to a pair of 
orthogonal axes 0x,0y with origin inside the area enclosed by Cp. 

For the Airy function F(x,y) the following formulas are well 
known: 


5 2 F _ d 2 F 

~ By 2 Y y c)x2 


a 2 F 

Xy = Y x - - 5x^y 


( 1 ) 


with the usual significance of the symbols. 


For any curve l the tangent and normal vectors denoted by t and 
n, respectively, are oriented so that the pair t,n is superposable in 
a rigid displacement of the plane on the oriented pair, -x,y. 

There is then obtained 




(2) 


If X n , Y n denote the components of the specific force corresponding 
to any point of l and in the oriented direction of n, 4 from the for- 
mulas 


Xn ~ X x n x 
y n = X y n x 


+ X y n y 
+ y y h y 


(3) 


Specifically, X n , Yn refer to the actions that the molecules situ- 
ated on the issuing side of the boundary exert upon the molecules on the 
other side. 
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and from equations (l) and (2) there is Imm ediately obtained 



If M n denotes the moment with respect to the third axis of 
coordinates z, of the vector of components X n ,Y n , there follows from 
equations (l) to (3) 


M n = x Y n - y X n = 


F - 


BF 
X 51 


- y 


Bf 

Sy 


(5) 


In the following discussion, c will denote any closed curve con- 
tained in T and inclosing Cq. When n is identified with the outer 
normal, the positive direction of transversing the curve, that is, the 
direction of t, is then such as to leave the area to the left. 

By making l coincide with c there is obtained from equations (4) 
and ( 5 ) 5 




x 


Bf Bf 
Bx “ y By 


M 


(7) 


^Formulas (6) and (7) are in accordance with a property of the Airy 
function already brought to light. (See reference 5.) 


1 1 ii ill i im i ii i na i mu 
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where X and Y denote the components of the resultant of the forces 
transmitted across the curve c from the interior toward the exterior 
and M is their resultant moment with respect to z. 

In the absence of body forces, however, X,Y,M are the corres- 
ponding components of the resultant force and moment of the external 
stress acting on the boundary Cp of the plane system. It is clearly 
shown in equations (6) and (7) that, in the presence of external stress 
of a quite general type, the expressions contained under the differ- 
ential sign are necessarily many-valued. It follows that the function 
F(x,y) and the first derivatives will in general be many- valued, but 
second derivatives will be single-valued (on account of their signi- 
ficance as stresses). 

If the two functions satisfying equations (6) and (7) are denoted 
by i|r^, function \|rp-i| r 2 > together with its derivatives, is 

necessarily single-valued, which is sufficient to assure that F is 
of the type 


F(x,y) = cp(x,y) + ilr'(x,y) (8) 

with cp (x,y) single-valued and i|r’ (x,y) uniquely determined from 
equations (6) and (7). 

When 9 = arctan y/x, the many-valued function -M0/2n (harmonic 
and with single-valued derivatives) satisfies equation (7). It is then 
found immediately that \|r' (x,y) corresponds to the product of 9 by 
a linear function of x,y, the coefficients of which are determined by 
taking into account equations (6) and (7). 

There is thus found 

+'(x,y) = - ^ 9 (9) 

where 

M p = M + yX - xY (10) 

is the resultant moment with respect to the pole P(x,y) of the 
external forces applied to Cp. 
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On the basis of equations (8) and (9), the Airy function can now be 
given the form 


M r 


F(x,y) = cp(x,y) - rj Q 


(H) 


2. Decomposition of the Single-Valued Part of the Airy Function 

It will be very useful in the following discussion to adopt for the 
function cp(x,y) the decomposition formula of Poincar^ (see paper by 
Fichera already referred to) and to write 


cp(x,y) =cp Q (x,y) + cp^(x,y) + 



+ |3x + ry + S 


3 


log p + 


a cos 20 + b sin 20 (12) 

The functions cp 0 (x,y), cp^x^y) are biharmonic, the first regular 
at all points not external to the region bounded by Cq, and the second 
regular at all points not internal to the region bounded by Cq and 
converging at infinity. 

In agreement with the notation adopted in reference 3 referred to. 


<5 L.u,v>c) 




ds + 



D 


. du 

Av Au 

dn 


c 



(13) 


where, if u,v are biharmonic in T, the value of $£u,v;c] does not 

depend on the curve c, and the constants a, (3, y, 5, a, and b are 
expressed by the formulas 


a =$[cp,l; c] ^ 


8 = $ [cp, p 2 ;cj > 

P = - 2<j>[cp,x;cj 

) 

a = cos 20; c] \ 

T = - 2 $ [cp,yJc] ^ 


b = |-$[cp,p2 sin 20; c] ^ 


( 14 ) 
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3. Expressions of Components of Displacement 

If e x , ey, e X y_ are the strains, cd the rotation, u,v the com- 
ponents of the displacement, following well-known formulas^ may he 
written: 


cSu 


5u 

l 

5cd 

1 ^ e xy 

.50 

dx 

e x 

5y ~ 

- eyrj - CD 

dx 

2 5x 

Sy 

> 

bv _ 

1 

dv 


ScD 

Sey ^ 2. 

S e xy 

dx 

— e-x-v + cd 
2 ^ 

Sy “ 

e y 

Sy 

dx 2 

Sy J 


(15) 

where e x , ey, e^y are connected with X x , Yy, Xy hy the relations' 7 


where 


^ - | (X x - T Y y ) 
ey * | (*y - T X x) > 
e xy = 2 ^ % 


(16) 


1 

in 

the 

case 

of plane 

stress 


l-p2 

in 

the 

case 

of plane 

strain 

(17) 

u 

in 

the 

case 

of plane 

stress 


v/ (l-v 2 ) 

in 

the 

case 

of plane 

strain 

(18) 



6 (See, for example, reference 6.) 
^ (See, for example, reference 7.) 
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In equations (16) to (18), V denotes Poisson's ratio and E, 
Young ' s modulus . 

From equations (l), (15), and (16) there follows immediately 

du 1+V b*F 

5y E Sx§y 03 

> 

bv k /S 2 F 
Sy. E ^ x 2 " by 2)^ 


_ k 5 af 
5x E c)y 

5a> _ k ciAF 
by E Sx 

By denoting now as l any line internal to T that unites the 
points P Q = (x 0 ,y 0 ) and P = (x,y), the following expression may he 
written: 


(19) 


( 20 ) 


and 


Su k /3 2 F a 2 F) 

Sx ~ E U7 - Bx 2 11 


bv 

bx 


1+1? a 2 F 

e a x a y 


+ Oi 


aao 



CDdq = TJCO - q 0 tt > 0 



Tjdkjo 


( 21 ) 


where rj denotes either of the variables x,y. 

Bearing in mind by means of equations (19) to (21), there is 
obtained with some transformations 
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u(x,y) = u(x 0 ,y 0 ) - [y co(x,y) - y Q a>(x 0 ,y 0 j] + 

i f/[^-»i>-^> / hi)} 

v(x,y) = v(x 0 ,y 0 ) + [x <o(x,y) - x G ffi(x 0 ,y 0 )] + > (22) 
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4. Airy Function Corresponding to Plane Dislocation in 
Doubly Connected Field 

With reference to case (a) of the introduction, the system is now 
assumed subjected to a state of strain due to a plane dislocation. In 
the decomposition of the Airy function defined by equations (10) to (12), 
the constants a, 0, and y will be shown to represent the three 
characteristics of the plane dislocation as determined by the latter 
equation. 

In order to show this relation, it is first observed that from 
equation (22) the necessary and sufficient conditions for the single 
valuedness of the components of the displacement and of the rotation 
are clearly expressed by 
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PdA F , 

J dr as = 0 

c 



The functions cp Q (x,y) and^Pj (x,y) from their nature® satisfy equa- 
tion (23) . 

Because the functions of log p and 0 are also harmonic and 
moreover Because log p is single valued with its derivatives and 0 
is many valued, hut with first derivatives single -valued, the two terms 
behave similarly. 

It is not difficult to show, by substituting, that also the func- 
tions cos 20 and sin 20 satisfy equations (23). 


®This is clear for equation (23.1) on the basis of a known theorem, 
if it is assumed that the functions ACPq and AC|>]_ are harmonic. As 
regards equation (23.2), it is noted that on the basis of equations (2), 

/ dAF dy\ , SlAF , / 3AF a _\ , 

a y 'sr ^ - ( y w - a ) **■ 

from which, by taking into account the biharmonic character of F the 
expression on the left of this equation is proved an exact differential. 
The left member of equation (23.2) is of the form of integrals of exact 
differentials. From the fact that <Pq has for its field of regularity 
the region enclosed by Co> it is then concluded that the left side of 
equation (23.2) vanishes for F= cpQ. The same result is obtained for 
F = cp^_ if it is considered that its field of regularity is composed of 
all the points of the plane not internal to the region enclosed by C-^. 
Similar considerations apply to equation (23.3). 
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Because of the properties of the terms that enter in the expres- 
sion of the Airy function (see equations (10) to (l2))> only those 
contained in the expression 


F*(x,y) = (a p 2 + px + yy) log p + ^ [x Y - y x] 0 (24) 

give rise to many-valued displacements. 

The increments received hy u(x,y), v(x,y), and <n(x,y) are 

denoted hy u*, v*, and a)*, respectively. When a complete circuit 
(started from the point of coordinates x,y) is effected about the 
hole, equations (22) that Eu>*/k, E(u* + ya>*)/k, -E(v* - xm*)/k 

agree with the values assumed hy the first members of equation (23) 
for F(x,y) = F*(x,y). 

If this result is taken into account, it is found from equa- 
tions (23) and (24) that 

U* - - 2» | [4 ay + 2r - y] 

V* . 2 ! I [4 (K t 2 P - x] 


* ok 

a> = 8< - a 
Jl 

Because a state of strain with surface forces (except for body 
forces) at the boundaries corresponds to each dislocation, there must 
be set X = Y = 0. From equations (25) the result is then obtained; 
The three characteristics of the plane dislocation (constants of the 
cut), are given by 




(25) 


J 


- l r 

- 8 *| a 
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The first two characteristics correspond to parallel fissures and the 
third to a radial fissure^. 

These fissures are assumed to be assigned on the basis of the par- 
ticular dislocation considered so that, when they are denoted by 2, 
m, and r, the following expression is obtained: Corresponding to 

a general plane dislocation of characteristics l, m, and r the 
Airy function in a doubly connected field is capable of the decompo- 
sition 

F = cp 0 + cp x + JL £ p 2 + mx - 2yl log p + 


6 log p + a cos 20 + b sin 20 


(26) 


The unknown part F = F - F* of the right-hand side of equa- 
tion (26) is uniquely determined by the condition that the forces on 
the complete boundary of the plane system are zero. 

Explicitly on the basis of equation (4) this condition is given 
by the equations 

> on C Q and Cj_ 

J 


bF 


5f* 

dx 


+ const. 


5f 

by 


bF* 

bj 


+ const. 


Q 

It is evident that these definitions give an immediate generaliza- 
tion of those used in the work of V. Volterra in the case of a field 
having the form of a circular ring. That is, the dislocation will be 
said to correspond to a uniform fissure when the displacement that 
generates it is a translation and to a radial fissure when the displace- 
ment is a rotation. 
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5. Airy Function in the Case of Single -Valued Displacements 
in Doubly Connected Field 


All the elements are now at hand for determining the structure of 
the Airy function in the absence of dislocations corresponding to any 
external stress. The coefficients a, 3> and T will now also be 
determined. In fact, the conditions of single -valuedness in equa- 
tion (23), which are more synthetically expressed by 

u* = 0, v* = 0, (jQ* = 0 


on the basis of equation (25) lead to 

oc = 0, r = (l-T)Y/4jt, 3 = (l-T)x/4jt 

From these relations and from equations (10) to (12), the following 
result is obtained: The Airy function in the case of single -valued 

displacements in a doubly connected field is capable of the decomposi- 
tion 


F = cPq + cp-j_ + 6 log p + a cos 20 + b sin 20 + 


_ 1 _ 

2jt 


(X x + 


Y y) 


- ~ T - log P - M p 



(27) 


Remark. It is sufficient to glance at equation (27) in order to 
derive the necessary and sufficient condition that the Airy function 
of a doubly connected system free from restraints does not depend, in 
the absence of dislocations, oh the elastic constants of the system 
and that the external forces applied at the boundary of the hole have 
a zero resultant. The same is naturally true for the stresses. 


6. Structure of the Airy Function in Multiply Connected Systems 

The preceding considerations will now be extended to the case of 
a plane system of any number of connections. The doubly connected 
system already considered will be referred to a pair of parallel axes 
coinciding with those already adopted but having their origin at a 
general point of coordinates -|, -T), which may be external to the 
whole. It is evident that the previously derived formulas, if referred 
to moving axes (denoted by x,y) , are written by putting in place of 
x,y, respectively, x-g, y-Tj. 
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In particular the Airy function will have the expression (see 
equations (10) to (12)) 


F = cp Q + 



(x-£) + r* (y-n) 



+ a* cos 20 + 


b* sin 20 - i Mn 0* 

2jt 


(28) 


where p* is the distance between the points of coordinate x,y and 


£ , t] and 0 

coordinates 
6*, a*, and b 
x-|, y-T] are 


is arc tan zia. It should be remarked that the point of 


* 


x ~i 

,T] is interior to the hole. The constants a*, p *, 
are defined by the same equations (14) except that 
substituted for x and y and p* for p. 


> 


Equation (28) gives a decomposition of the Airy function referred 
to a pair of axes x,y completely arbitrary but with having the 

significance of coordinates of a point Q within the hole; this equa- 
tion is generalized to the case of plane systems with any order of 
c onne c t i vity . 


The region of definition of the Airy function will now be con- 
sidered as a region T bounded by the curve Cq and having n holes 
bounded by the curves Cp, 0 %, . . . C n ; |p, rjp (i = 1, 2, . . . n) 
will be taken as the coordinates of a point within the hole 

bounded by Cp (i = 1,2, . . . n) . 

Setting 

p i 2 = ( x - U 2 + 


0 - = arctan ^ ^ 

x -tl 

the Airy function in this region, by evident generalization of equa- 
tion (28), has the following form 



16 


NACA TM 1290 


n n 


F - *0 + £ + £ 


[“! p\ + Pi (x-Ej) + Ti + PJ log Pi + 


i=l 1=1 


i. cos 26 . + Id. sin 29 . 1 - ~ Y' M 0, 

i i 1 1 [ 2 it /_/ P 1 


(29) 


i=l 


In this equation, q>Q is biharmonic and is regular in the manifold 
bounded by Cq} <Pp (i = 1,2, . . . n) is biharmonic outside the hole 
bounded by Cp (i = 1,2, . . . n) and converges at infinity. 

If Xp, Yp ( 1 = 1,2, . . . n) are the components of the result- 
ant of the external forces applied at the boundary of the hole bounded 
by Cp and Mp is their resultant moment with respect to the axis 

oriented as z through Qp, = (y-Tlp)Xp - (x-|p)Yp + Mp denotes 

the resultant moment applied at Cp with respect to the axis oriented 
as z through the point P = (x,y). The constants ap, pp, yp, 6p, 
ap, and bp are determined by the formulas 


a i 

= ® 

[cp 

> c i] 


Pi 

= - 

2® x-|p, 

=3 

r± 

= - 

2 * y-tipd 

°D 

Si 

1 

3 

P 2 i> c i) 


a i 

1 

2 

® 

[ 3 , p 2 ± cos 

20p; cj 

Si 

H|<U 

II 

® 

[), p 2 ± sin 

20 i' op] 


( 30 ) 


where cp denotes the single-valued part of the Airy function. 
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It will new be shown that in equation (29) the coefficient^ of 
log does not depend on t^, (i = 1,2, . . . n) . 

Because of the linearity of the operator $, the following expres- 
sion may be written, for any value of i 

p ± = - 2 4> Qp, x; c-f] + 2 { ± $ [cp, 1} c[] = 0^ + 2 a ± % ± (31) 

where 

Pi = - 2 4> [cp, x; c j 

Similarly 

*i = f i + 2 “l \ 

with 

fi = - 2 $ [q3, yj c£] (34) 

Moreover 

S ± = ® [qp, P 2 ; cj + (li 2 + rii 2 )® Qp, 1; cj - 

2 li* [qp, x; cj - 2 •n ± <!> [cp, y ±} c[] (35) 

where P denotes the distance from the origin of any point P = (x,y). 

From equations (30.1), (32), (34), and (35), the following expres- 
sion is immediately obtained 

5 i = &i + “i (Ip 2 + n± 2 ) + F ± li + 7 ^i (36) 


■*-9a similar observation can evidently be made with regard to the 
coefficient of log p* In equation (28) . 


(32) 

(33) 
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with 

8 = $> [q>, P 2 ; cj (37) 

On the basis of equations (31), (33), and (36), the following equa- 
tion can he written: 


H Pi 2 + ^i ( x ~£i) + r± (y-n ± ) + 5 ± = a* p 2 + x + y + 5 ± (38) 

Equation (29) can therefore be written as 

n n r~r- ^ . 

F = cp Q + ^ cp± + ^ ) [^i p + Fi x + T i y + 6jl log P A + a i cos 26.^ + 
i=l i=l (■ 

b ± sin 20i - ^M p (l) 0i| (39) 

The meaning of cc^, p^, T ^ remains to be determined. For this 
purpose, u s *, v s *, cog* (s = 1,2, . . . , n) denote the increments 
received by u, v, <x> when a complete turn is made about the hole 
bounded by C s (s = 1, 2, . . . , n) starting from the point P s (x,y). 
It is evident that, of the terms of the expression equation (29) of the 
Airy function, only those of index s contained in the various sums 
can make a nonzero contribution to the formation u s *, v s *, o) s * 

(s = 1,2, . . . , n) , because the others are regular with all the deriv- 
atives in the region bounded by C s (s = 1,2, . . . , n) . 

Substantially, the same relation holds ture as in a double con- 
nected system and the expressions of u s *, v s *, o> s * are obtained from 
those of u*, v*, u>* (equation (25)) by simply substituting a s , p s , 

T s , X s , Y s for a, p, y, X, Y, respectively, and x-| s , y-Ti s for 
x,y. 
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u s * - - 2n | ^4 Og (y-T] s ) + 2 V s - Ysj 

v s = 2it g- ^4 a g (x-| B ) + 2 p s - xJJ, (s = 1,2, • • • > n )^ 


“s’" = 8* | «b 


(40) 


It is then found that u B *, v s *, ao s * do not depend on the points 

(£s> %)• 

In fact, on the basis of equations (31) and (33), equations (40) 
may be written 


Ug = - 2jt — |-* v*g 


I E 


y + 2 r s - 


2 It 


1y 3’ 


'N 


v s* = 2jt | [1 a s x + 2 p s - xJJ, (s = 1,2, . . . , n) l (41) 


oj s * = 8jt — a E 
E 


j 


from which immediately follows the result: If 1 m-p and r^ are 
the characteristics of a plane dislocation corresponding to a cut that 
goes from C 0 to Cj_ (i = 1,2, . . . , n) , 


Ctj[ = - -JL v± 

^ 8jtk 1 




Pi = 


E 

4jtk 




(42) 


T. = _ JL l • 

1 4rtk 1 


J 


It is also clear that the 0^, Yi correspond to parallel fissures 
and the to radial fissures. 
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From equations (38) and (42) the final result is thus obtained: 
The structure of the Airy function corresponding to a general plane 
dislocation in an (n+l)-fold connected field is defined by the 
formula: 


F = cp Q + 


n 

S>± 



r i —2 

X P + ULs X 
2 ^ L 



+ 



n 

log Pi + 52 
i=l 



cos 


2 + bp 



In order that the displacements be single -valued, it is evidently 
necessary and sufficient that 



Vi* - 0 
cju^* = 0 


A 

> (i = 1,2, 

J 


n) 


From equation (41) : The necessary and sufficient condition for 
the single-valuedness of the displacement in an (n+l)-fold connected 
plane system is that the following equalities be satisfied: 


ci^ = 0 




r = TF Yi 




(i = 1,2, 


J 


n) 


and therefore the Airy function in the absence of dislocations in an 
(n+l)-fold connected field is capable of the decomposition-^- 


-^From (43) there clearly results the condition (analogous to that 
expressed in the remark of section 5) of the independence of F on the 
elastic constants of the system. 



NACA TM 1290 


21 


sin 20q] + 


|(X ± x + Yi y) log Pi - Mp^ 0i| (43) 

Remark. Equation (42) evidently also gives the expression of the 
Airy function in the presence of forces or moments concentrated at 
interior points of the plane system. If, for example, at the point 
Q = (|*,tj*) there acts a concentrated stress (concentrated force or 
moment or both systems), one of the terms of the stun in equations (42) 
corresponds to it. If this term of the sum is, for example, of index 1, 
<Pl is biharmonic at the exterior of each circle of center Q, and cp-^ 
converges at infinity, then 8-^, aq, bq are always defined by equa- 
tions (37), (30.5), and (30.6), respectively. In these relations, cq 
is any circumference of center Q interior to Cq and such as to leave 
the regions bounded by Cq, Cg, ... on the outside; Xq, Yq, M p (1) 
denote the components of the resultant of the concentrated stress at Q 
and the resultant moment with respect to the straight line oriented as 
z through P = (x,y). 


f = <? 0 + S ^i + S . [ 5 i 1o s Pi + a i cos 2 0q + bq 


7. Mechanical Significance of Coefficients 8, a, b 

As evident from the heading of the present section, the case of a 
doubly connected field is considered for simplicity, but the extension 
to the case of any number of connections will be evident. 

The coefficients 8, a, b are not determined, as are a, (3, 
and Y, by the condition of the problem; they nevertheless have a 
mechanical significance that is now given. 

It will be shown that these coefficients are expressed in terms of 
characteristic elements of the deformation and of the stress such as 
the rotation, the coefficient of surface dilatation, and (only for the 
coefficient 8, however) the astatic coordinates of the external stress 
acting on one of the boundaries of the plane system. 


22 


MCA m 1290 


It is first necessary to establish several preliminary formulas. 
We set 


F** = ^[?Y-yX-Ml]e = -||e (44) 





5** = 


'•*, p 2 > < 

a 






a** = 

I* 

&**, P 2 

cos 

20} 

c] 

► 

(45) 


b** = 


[f**, p 2 

sin 

20} 

<3 , 



From equations (44) 

and (45) 

, if 

equation 

(13) 

is 

taken into 

account. 

there is obtained 










5** 

1 

" 2* 

§o Its o l ] 




(46) 



a** = 

= b** = 0 





(47) 


where Xg, yQ denote the coordinates of the point Pg from which the 
path on c begins. 

On the basis of equations (10), (ll) , (14) , (44), and (45), and 
when the linearity of the operator $ is taken into account, the fol- 
lowing relation may be written: 


6 = ft [f, P 2 J c] - 5** 

where as when equation (47) is also accounted for, 

a = — $ [F, P 2 cos 20; c] 


(48) 




( 49 ) 


b = — 0 |F, P 2 sin 20} c| 
2 
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Taking into account equations (2), 


c c c 

/C xd l- ya l) (so) 

c 

From equation (50), on the basis of equations (4) and (6) there is 
immediately obtained 

(x^ + y Y n ) ds (51) 

c c 

When two of the astatic coordinates (which two are required is 
evident) of the stress acting across c are denoted by a-^ and a-yy> 

equation (51) can be (see reference 8) written (see also equation (46)) 


ds = x Q X + y 0 


Y + 


f da = 2 k 5** + a xx + 


®yy 


(52) 


From equation (22.3) there is immediately derived 


E . d£F , 

r- ds = - 5 — ds 

k dn 


(53) 


and also 


E 




(54) 
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where D denotes the divergence of the displacement u,v (coefficient 
of surface dilatation) . 

Mechanical significance of 8 and of its invariance with respect 
to c. - When equation (48) is expressed explicitly, and on the basis of 
equation (13), the following expression may be written: 


6 " j J p2 TGT ds + 4 J I ds " IT ds ' " 

Taking into account equations (52) to (54), 


& = 


E 

8jtk 


P 2 dn> - _iL_ / D ds > + -i_ (a xx + a-,™-) 

1-T J dn [ 2jt yy/ 

c 


(55) 


or, taking into account equations (52) and (54) 


8 = 


E 


P 2 ® ds - 


8rt k (1 -t) J dn 


J D TET dS | + h ^ axx + 


Equations (55) and (56) give the required expressions of & as a 
function of the characteristic elements of the strain and the stress^. 

A closed curve c' of the same type as c is now considered and 
the area of the band bounded by c and c' is denoted by A cc '. 

Also denoted by Oyy are the two astatic coordinates of the 

stress that acts on the band across the curves c and c* and is 
defined by 


-'- 2 If the state of strain is due to a simple plane dislocation on one 
of the boundaries, then a^x = ciyy = 0 and it is sufficient to eval- 
uate 8 corresponding to or ctyy to obtain a simpler expression 

than equations (55) or (56) . 
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a xx = / x X n ds - (x X n ds 


ayy = 1 7 ? n ds “ 7 *n ds 


*\ 


( 57 ) 


J 


When 8‘ denotes the expression of 5 relative to c', there is 
obtained from equations (56) and (57) 


8 - 8 * = 


E 


8 it k (!■ 




S “x + °y) ds ' 2 J D ( x *x + y “y) ds 


' 2 (J^ nx + §y dS + 2 J D ( x n x + y °y) ds > + ^ (°xx + Oyy) 


p V = n x + Sr 


(58) 


When the transformation formula of line integrals into surface 
integrals is applied to equation (58) and the harmonic character of D 
is considered 


D d A cc - + — (oxx + c^) ( 59 ) 

If I denotes the mean value of the linear invariant of the 
stress I and the independence of 8 on c in A cc i is accounted 
for, there is derived from equations (54) and (59) 


8 - 8 ' = 


E 


2jt k (t-1) 


cc 


I _ a xx + a yy 
A c c' 

that is, the invariance of 8 with respect to c expresses the prop- 
erty of the mean stress that assigns the mean value of the linear 
invariant of the stress^ 3 . 


dd For the properties of the mean of the stress mentioned above see 
reference 9. 
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Significance of a, b and of their Invariance with respect to c. 
On the basis of equations (13) and (49), ’ — — * - 


a = 


16rt 


r\ 


, 2 2\ dAF - 

(x - y > -as ds 




i / 


y (so) 


b = -J AF ^ ds 

[ c c 


When equations (53) and (54) are considered, equations (60) reduce to 


a = 


E 


16itk 


U 2 - y 2 ) to - j4t 


h = 


ffiik ^ J xy do3 - ds 

c c 


> ( 61 ) 


or 


a = 


h = 


E 


2 , r 2\ ® ^ /d(x 2 - y 2 ) 


i a s(TT T r ) J (x - y } aE ds - 

c 


dn 


D ds > 


E 


8itk(l - T) 


d D , 

^ dE ds 


dxy 

dn 


D ds 


y ( 62 ) 


With the aid of equations (61) and (62), a and h are thus 
expressed as functions of the rotation and of the coefficient of sur- 
face dilatation. It is also found (from equations (60)) that the 
invariance of a,b "with respect to c is derived from a known prop- 
erty of harmonic functions. 

On the basis of equation (l), the following expression may be 
written: 


t(x ' y) ds = / t(x ' y) |_ “x + \-£ - ~sk) “yj 


/s y ax. 


ax x ax v \ 


ds 


( 63 ) 
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where i|r(x,y) denotes any integrahle function of x,y. From equa- 
tion (63) there is obtained 




a cS Xy 

5y 


5 IjF Xy 

5x 


+ 



CC ' 



, d 2 i if 
y dx 2 


-2X y 


d 2 ^ 

dx dy 


d A, 


cc ' 


+ 



a z x y 

dxdy 


d A, 


cc 1 


(64) 


Taking into account that equation of congruence of Saint -Venant 
which remains in effect in the case of plane systems expressed in 
terras of the components of the stress it is seen that the last integral 
of the second member of equation (64) is zero. 


From equation (64) there is then obtained 


I , dAF . . dAF , 

^ ds *,ris ds = 


St Y Si|r \ 

Y y " 5y x yJ ** 


J| x * “ al x yJ \\ 


ds - 


5 


St Y St „ ) 

32 Y y ■ 32 x y/ x 


(|f X x - 32 0 n y_ ds " 


* 


x + y 2x 

ay 2 y ax 2 y ' 


^ ^CC* 


( 65 ) 


By taking into account equations (l) and (3), the following expression 
is obtained from equations (60) and (65) for t = x 2 - y 2 
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- a ’ = 55 { f( 7 *n ~ xX n) ds -Jlfl -. n - xXj ds - f (Y y - X*) d A cc , 
l c e« A.., 


( 66 ) 

where evidently a 1 denotes the expression of a in correspondence with 

c* . 

By taking into account equation (57) and denoting by X 2 , Y y the 
mean values of X^., Y y in A cc t, the following expression may be 
derived from the invariance of a and from equation (66): 


Z X - Yy = 


a xx ~ a yy 

K 


cc ' 


from which it is seen that the invariance of a with respect to c 
expresses one of the properties of the mean stress. 

From equations (60) and (65), for tJ/ = xy there is obtained, on 
the basis of equations (l) and (3) 


b - b‘ 



-/'(* Y n + yXj ds + J\x Y n + yXj ds + 2 / ly d A cc , 
c c * A cc' 


where the meaning of b* is clear. 


(67) 


The astatic coordinate of diverse indices of the stress acting on 
the band bounded by c,c* across its boundaries is denoted by a^y 
and the mean value of Xy in A cc r is denoted by X y . From the invar- 
iance of b and from equation (67) there is immediately obtained, if 
the second fundamental equations of statics are taken into account. 


a. 




It is therefore seen that the invariance of b with respect to c 
expresses that property of the mean stress which assigns the mean value 
Of Xy. 


1340 
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Eemark , If c is assumed to coincide with a circumference^ of 
center 0 and of radius E, equation (55) becomes 

2rt 

8 - He[ “* - d9 ] + h + V> <68) 

and presents 6 as a linear combination of the increment of the rota- 
tion, because of a turn about the hole, of the mean value D of D on 
the circumference c and of two astatic coordinates. In particular if 
(0* = 0 equation (68) becomes 

•^2 ^ 

6 = " T k(l-T) ~ a xx " a yy (69) 

and: In the absence of dislocations and in the presence of disloca- 
tions that permit single-valued rotation (a = 0), 6 is equal to a 

linear combination of the mean value of the coefficient of surface 
dilatation along each circumference of center 0 and of two of the 
astatic coordinates of the stress that is transmitted across the cir- 
cumference. This linear combinatidn is therefore invariant with respect 
to the circumference. 

In particular for a plane system with circular hole subject to a 
dislocation that leaves the rotation single-valued (a = 0), it is 
sufficient from the fact that c is the circumference which bounds the 
hole and that correspondingly a YY = a^y =0 to derive the result that 
5 expresses, except for an obvious coefficient, the mean value on c 
of the surface dilatation. 

The same relation holds in the absence of dislocation, provided 
the external stress applied at the boundary of the hold satisfies the 
condition a yy = =0 (in particular if the hole is free from exter- 

nal stress). 


•^Naturally for the validity of the remark it is required that 
there exist a circumference entirely internal to T . If this is not the 
case, the biharmonic function F does not lose significance in a region 
T* comprising T in which circumferences can be drawn; the remark 
retains its validity provided the elastic system is considered extended 
in T* . 
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8. General Expressions of Components of Displacement 
in Doubly Connected Field. 

From the preceding developments, it is easy to obtain the expres- 
sions for the components u,v of the displacement. There is thus the 
possibility of considering problems in which the displacements are 
assigned on the boundary or mixed problems. 

For simplicity, the case of a doubly connected system will be con- 
sidered, but it is not difficult to extend the results obtained to sys- 
tems of any order of connectivity. 

The values assumed by u,v are denoted by u,v when, in equa- 
tions (26) and (27), CPq/P-^ are assumed to be zero in addition to 5, 

a, and b. The expressions of u,v are obtained from equation (22) by 
substituting in them F - cpQ - cp^ - 8 log p - a cos 20 - b sin 20 in 
place of F. When a rigid displacement is singled out and expres- 
sion (26) is assumed for F itself, there is obtained in case (a): 


= 2 £0 + iil sin 20J + m|~ill log p - Jktl cos 20 + — llj 


rp J^0 sin 0 + cos 0 log p + i cos ©j | 

= ~{ 1°S P “ "X 1 cos 29 + “ x 1 sin 20J 

[ T-l T+l ”11 

- 0 cos 0 + — — sin 0 log p + —j— sin 0 ^ 


and in case (b) when the expression (27) is assumed for F, 
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u = 


k 

8«E 


[x [2 (t 2 -2T-3) log p + (1+T) 2 cos 20 + (l-T) 2 - 2(3 +t)] + 


"N 


Y(l+T ) 2 sin 20 - 4(1+T) M 


sin 0 


) 


v = 


8«E 


|x(l+T) 2 sin 20 + Y |2 (t 2 -2T-3) log p - (1 +t) 2 cos 20 + 
(l-T) 2 - 2(3+T)] + 4(1+ T) M ££ ^ J 


J 


Next, if u*,v' denote the values assumed hy u,v when F is 
identified with 5 log p + a cos 20 + h sin 20, there is obtained on 
the basis of equation (22) 


u' = ^-8(1 +t)cos 0+4a cob 0(cob 2 0- t Bin 2 B)+Zb [(3+t)cos 2 0+(l-T)sin 2 0]ain 0 

v' = |L |-8(l+f)sin 0+4a Bin 0 (Tcob 2 0-Bin 2 0)+2b [ (3+T)sin 2 0+(1 -t)oob 2 0]cob 0 


Finally, if u ,, ,' v ■ ,, denote the values assumed by u,v when for F 
there is assumed cPq + 

u = u + u’ +u’ 1 
v = v + v’ + v 1 1 

where u’^v* ’ are given by equation (22) for F = 

It may be useful in many problems to assume for <Pq and a 

special series expansion in trigonometric binomials. By known theorems, 
the biharmonic functions ^Pq, can be developed into uniformly con- 

verging series, the first at each point at infinity, and the second for 

p > 0 



hIo. i-ila 
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where 


T . = a. cos n 0 + b. sin n 0 

in in in 


(i = l,2,3,4;n = 0,1,...) 


(70) 


(71) 


In the series (70.2), the term T 42 of the type a cos 20+b sin 20, 

already separately considered, has been suppressed. It is sufficient 
to substitute equation (70) into equation (22) to obtain 


n= 2 , 3 , . . . 


u" = u" + 
20 


v" = v" + 
20 


E u "n + E* u " 


in 


i=l 

3 


i=l, . . .4 


n= 2 , 3 , . . . 
* 


L • l 

i=l i=l, . . .4 


v". 

in 








(72) 


with 
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u "ln = - I ( T +l)p T1 - 1 [^ T ln cos 0 - Bln 0 


(n = 1,2,...) 


e 

r 


"in = “ I ( T+1 )p n_1 £ T in sin 0 + cos 0 J 

^ [j + ( T+1 )] sln 9 ‘ \ T | n “ &-(°.+2)(T+l)3 T 2n cos 0 

" 2 n = “ 1 P n+1 | [| + (T+l)]cos 0 - [Mn + 2)(T+1)] T 2n sin 0 


u "2n - 5 P n+1 


(n= 0,1,2, . ..) 
k / T .ii.-(n+l)J^ n T 3ri cos 0 + sin 0 


1(73) 


u "3n ~ ^ ( t+1 )p 


a T 3n 

d 0 


(n — 1,2,...) 


‘] 

v" 3n = | (T+l)p-( n+1 )[^ T 3n sln 0 - cos 0 

s' 

u "4n = 1 P 1_n | | + (T+l)Jsin 0 +jl-(2-n) (T+lfjcos 0 T 4n 

v "4n = | P 1_n | [n ' ( T +l)Jc°s 9 -^g 11 + &- (2-n) (T+l)J sin 0 T 4n 




(n — 3,4, . . . ) 




The asterisk on the summation signs denote that the term of index 4,2 
is not considered. It is^necessary moreover to consider that for n = 0 
the terms that contain ~2Q are suppressed in equation (73). 


Translated by S. Reiss 
National Advisory Committee 
for Aeronautics. 
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